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The exact solution to the Schrodinger equation for the rigid body with the 
given angular momentum and parity is obtained. Since the quantum rigid 
body can be thought of as the simplest quantum three-body problem where 
the internal motion is frozen, this calculation method is a good starting point 
for solving the quantum three-body problems. 
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1. INTRODUCTION 

The three-body problem is a fundamental problem in quantum mechanics, 
which has not been well solved. The Faddeev equations [1] provide a method 
for solving exactly the quantum three-body problems. However, only a few 
analytically solvable examples were found [2]. The accurate direct solution 
of the three-body Schrodinger equation with the separated center-of-mass 
motion has been sought based on different numerical methods, such as the 
finite difference [3], finite element [4], complex coordinate rotation [5], hy- 
perspherical coordinate [6-8], hyper spherical harmonic [9-11] methods, and a 
large number of works [12-16]. In those numerical methods, three rotational 
degrees of freedom are not separated completely from the internal ones. In 
this letter we present a method to separate completely the rotational degrees 
of freedom and apply it to the quantum rigid body as an example. 

The plan of this letter is organized as follows. In Sec. 2 we shall introduce 
our notations and briefly demonstrate how to separate the rotational degrees 
of freedom from the internal ones in a quantum three-body problem. The 
exact solution to the Schrodinger equation for the rigid body with the given 
angular momentum and parity is obtained in Sec. 3. A short conclusion is 
given in sec. 4. 

2. QUANTUM THREE-BODY PROBLEM 

Denote by rj and by Mj, j = 1, 2, 3, the position vectors and the masses of 
three particles in a three-body problem, respectively. The relative masses are 
nij = Mj/M, where M is the total mass, M = YJ Mj. The Laplace operator 
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in the three-body Schrodinger equation is proportional to J2j=i m J A Tj , 
where A rj is the Laplace operator with respect to the position vector rj. 
Introducing the Jacobi coordinate vectors x and y in the center-of-mass 
frame, 

x = -\ ; r i> y = \ ; (r2-r 3 ). (1) 

we obtain the Laplace operator and the total angular momentum operator 
L by a direct replacement of variables: 

3 

A = E m J lA rj = A x + A y , 

J'=l 

3 (2) 
L = 22 - ihr j x Vr, = L x + Ly, 

j'=i 

L x = -ihx. x Vx, L y = -ihy x y y - 

The three-body Schrodinger equation with the separated center-of-mass mo- 
tion becomes 

- (V/2MJ {A x + A y } * + VV = EV, (3) 

where V is a pair potential, depending only upon the distance of each pair 
of particles. 

In the hyperspherical harmonic method [11], for example, two Jacobi 
coordinate vectors are expressed in their spherical coordinate forms, 

x ~ (pcosu,0 x ,ip x ), y ~ (psmu,0y,ip y ). (4) 

where p is called the hyperradius and 8 x , tp x , y , <p y ) are the five hyper- 
angular variables. The wave function is presented as a sum of products of a 
hyperradial function and the hyperspherical harmonic function, 

There is huge degeneracy of the hyperspherical basis, and the matrix ele- 
ments of the potential have to be calculated between different hyperspherical 
harmonic states [10], because the interaction in the three-body problem is 
not hyperspherically symmetric. 

The quantum rigid body (top) can be thought of as the simplest quantum 
three-body problem where the internal motion is frozen. To solve exactly the 
Schrodinger equation for the rigid body is the first step for solving exactly 
the quantum three-body problems. Wigner first studied the exact solution 
for the quantum rigid body (see P.214 in [17]) from the group theory. He 
characterized the position of the rigid body by the three Euler angles a, (3, 
7 of the rotation which brings the rigid body from its normal position into 
the position in question, and obtained the exact solution for the quantum 
rigid body, which is nothing but the Wigner D-function. For the quantum 



2 



three-body problems, as in the helium atom, he separated three rotational 
degrees of freedom from three internal ones (see Eq. (19.18) in [17]): 

*<m(ri,r 2 ) = 2 D^(a,/3,7)V v (ri,r 2 ,w), (5) 

V 

where ri and r 2 are the coordinate vectors of two electrons, uj is their an- 
gle, and the Wigner D-function form [17] has been replaced with the usual 
D-function form [18]. Wigner did not write the three-body Schrodinger equa- 
tion explicitly. As a matter of fact, the three-body Schrodinger equation (3) 
becomes very complicated if one replaces two coordinates vectors of electrons 
with the Euler angles as well as n, r 2 , and oj for the internal motion. On 
the other hand, Wigner's idea, to separate the degrees of freedom completely 
from the internal ones, is helpful to simplify the calculation for the quan- 
tum three-body problem. Hsiang and Hsiang in their recent paper [19] also 
presented the similar idea. In this letter we will develop the idea of Wigner 
and obtain the exact solution of the Schrodinger equation for the rigid body 
without introducing the Euler angles directly. This calculation method is a 
good starting point for solving the quantum three-body problems [19,20]. 

The Schrodinger equation (3) is spherically symmetric so that its solution 
can be factorized into a product of an eigenfunction of the angular momen- 
tum L and a "radial" function, which only depends upon three variables, 
invariant in the rotation of the system: 

£i = x • x, £ 2 = y • y, & = x • y. (6) 

For the quantum rigid body, the potential makes the internal motion frozen 
so that those variables £j are constant. 

For a particle moving in a central field, the eigenfunction of the angular 
momentum is the spherical harmonic function Y^{d, tp). How to generalize 
the spherical harmonic function to the three-body problem without intro- 
ducing the Euler angles directly? As is well known, ^4( x ) = r ^^m(0>¥O> 
where (r, 9, <p) are the spherical coordinates for the position vector x, is a 
homogeneous polynomial of degree I with respect to the components of x, 
which does not contain r 2 = x x explicitly. 3^(x), called the harmonic 
polynomial in the literature, satisfies the Laplace equation as well as the 
eigen-equation for the angular momentum. 

In the three-body problem there are two Jacobi coordinate vectors x and 
y in the center-of-mass frame. We shall construct the eigenfunctions of the 
angular momentum as the homogeneous polynomials of degree £ with respect 
to the components of x and y, which do not contain £j explicitly. According 
to the theory of angular momentum [18], they are 

y) = E ^(x)34_V(y)<9, H,i-Q,m- i-q,L,m), 

M (') 

< q < £, when L = £, and 1 < q < i — 1, when L = I — 1. 

where (q, fj,, £ — q,m — fi\q, £ — q, L, m) are the Clebsch-Gordan coefficients. 
The remained combinations with the angular momentum L < £ — 1 contain 
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the factors £3 explicitly [20]. In other words, the eigenfunctions of the total 
angular momentum L 2 with the eigenvalue £(£+1), not containing the factors 
£j explicitly, are those homogeneous polynomials of degree £ or degree (£+1). 
Let us introduce a parameter A = or 1 to identify them: 

y ( L x)q ^y) = E^(x)^ +A (y) 

x {q, fi,£ — q + A, m — /j,\q, £ — q + X,£, m), (^) 
A = and 1, A < q < £. 

yfm^i^y) ls tne common eigenfunction of L 2 , L3, L 2 , L 2 , A x , A y , 
A xy , and the parity with the eigenvalues £(£ +1), m, q(q + 1), (•£ — q + 
A)(^ — q + A + 1), 0, 0, 0, and (— l) e+x , respectively, where L 2 and L3 are 
the total angular momentum operators, L 2 and L 2 are the "partial" angular 
momentum operators [see Eq. (2)], A x and A y are the Laplace operators 
respectively with respect to the Jacobi coordinate vectors x and y, and A xy 
is defined as 

A 92 ■ &2 &2 (9) 

xy dxidyi dx 2 dy 2 dx 3 dy 3 

Because of the conservation of the angular momentum and parity, the 
solution ^£ m A(x, y) of the Schrodinger equation (3) can be expanded in 

terms of J^^ - ^ 9 (x, y), where the conserved quantum numbers £, m and A 
are fixed. Since those equations are independent of m, we can calculate them 
by setting m = £, where [18] 



3#(*,y) 

[(2q + l)l(2£ - 2q + l)l] 1/2 




q\{£ - q)l2 e + 2 7r 



(2q + 1)\(2£ - 2q + 3)! 



(xi + ix 2 ) q (yi + iy 2 ) 



t-q 



(10) 



| 1/2 {( 9 -l)!(^-g)!2^ 2 7r} 1 



2q{£-q+ !){£+!) 
xi + ix 2 ) q ~ 1 (yi + iy 2 y~ q {(xi + ix 2 )y 3 - x 3 (yi + iy 2 )} X • 

By substituting \I/^(x, y) into Eq. (3), a system of the partial differen- 
tial equations for the coefficients can be obtained. The partial differen- 
tial equations will be simplified if one changes the normalization factor of 

34 m)9 (x,y), namely ^ +A)<z (x,y) in Eq. (11) is replaced by Q q x (x,y), 
which is proportional to ^ +A ^(x, y): 

e 

Wx,y) = £ (x,y), 

g=X 

Qf(x,y) = {{ q -\)\(£- q )\}-\ Xl+lX2 y-\ yi+l y 2 Y-i (11) 

X {(Xi + iX2)2/3 - ^3(2/1 + «/2)} A 

A = 0,1, A<g<^. 
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The partial differential equations for the functions Y>g A (£i, £2, £3) are: 
l A< + 4^ + 4(£-g + A) 1 ^ + 2( (? -A)^i 




(12) 



This system of the partial differential equations was first obtained by Hsiang 
and Hsiang [19]. It is a good starting point for solving the quantum three- 
body problems [19,20]. 

3. QUANTUM RIGID BODY 

For the quantum rigid body, the potential preserves the geometrical form of 
the rigid body fixed. It can be replaced by the constraints: 

£1 = const. £2 = const. £3 = const. (13) 

Therefore, the solution of the Schrodinger equation for the quantum rigid 
body can be expressed as 

*«A(x,y) = ]T /fQ?(x,y). (14) 

q=X 

where f^ x are constant. Recall that Qq A (x, y) is the solution of the Laplace 
equation. Due to the constraints (13) some differential terms with respect to 
£j in the Laplace equation should be removed so that the Laplace equation 
is violated, namely, the rigid body obtains an energy E. On the other 
hand, as a technique of calculation, we can calculate those differential terms 
first where £j are not constant, and then set the constraints (13). The 
contribution from those terms is nothing but the minus energy — E of the 
rigid body. 

In the calculation, we first separate the six Jacobi coordinates [see Eq. 
(4)] into three rotational coordinates and three internal coordinates. The 
lengths of x and y and their angle 00 are 

r x = r y = a/6, cosw = 

Obviously, those three variables are also constant in the constraints (13). 
Assume that in the normal position of the rigid body the Jacobi coordinate 
vector x is along the Z axis and y is located in the XZ plane with a positive 
X component. A rotation R(a, (3, 7) brings the rigid body from its normal 
position into the position in question. The Euler angles a, (3, and 7 describe 
the rotation of the rigid body. The definition for the Euler angles are different 
from that of Wigner (see Eq. (7) and Ref. [17]) because x and y here are 
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= cos a, 




= COS0J;, 




= sin a, 


Sec 


= sin 9 X , 



the Jacobi coordinate vectors. To shorten the notations, we define 

C/3 = COS P, Cj = cos 7, 

Cy=COS9y, C = cosuj, , . 

sp = sin (3, s 7 = sin 7, ^ ' 

s y = sinOy, S = sinuo. 

According to the definition, we have [18] 

SaCpC-f + C Q S 7 -S Q ,C /3 S 7 + C Q C 7 

-S/3C 7 SpS-, C/3 

x\ + ix 2 = r x e ta Sf3, yi + iy 2 = r y e ta (cpc^S + spC + is y S) , 
X3 = r x c/3, y 3 = r y (-spc^S + cpC) . 

Through the replacement of variables: 

{r x ,9 x ,<p x ,ry,6 y ,ipy) — ► (r x ,r y ,u,a,/3,^), 
a = fx, P = Ox, 

C = C X Cy + S X SyCOS((f X - <Py), (18) 
S X Cy C X Sy COS {<Px fy) 



(17) 



cot 7 



Sy Sm((p x - (fy) 



we obtain 



i d 2 

A x = — ~^ r x^ , 

r x or* 

A —— r | ——s— 

r y dr 2 y r 2 S duo duo 



(19) 



where the neglected terms are those differential terms only with respect to 
the rotational variables a, (3 and 7. Now, 




l —r+-—r 1 9 s 9 \ ip 

x r„ dr 2 y r^S duo duo \ ' 



r x v x ' y y 1 y 



(20) 



£*« A (x,y) 



£ } j=const. ' 



where *^(x, y) is given in Eq. (14). 

Through a direct calculation, we obtain 

f|U^V,y) = ^«?(x,y), 

l^ f(x , y) .(lzl±^±A±i) Q?(3£ , y) , 

rlSdu^du^ 1 ' yj (21) 
= {(* - 9) [(* - 9 + 2A ) cot2 w - !] 
+ A(cot 2 ^-l)}Qf(x,y)/r2 

- (g - A + 1) (2£ - 2q + 2A - 1) {C/S 2 ) Q^fry)/ (r x r„) 
+ (g - A + - A + 2)S~ 2 Qf + 2 (x, y)/ r 2. 
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Therefore, the coefficients f q satisfies a system of linear algebraic equa- 
tions with the equation number (£ — A + 1): 



(2ME/h 2 )f q x = {q(q + l)/r 2 x + (£ - q + \)(£ - q + \ + l)/r 2 y 

+ [{t-q){t-q + 2A) cot 2 CO - {£ - q) + A (cot 2 - 1)] /r 2 } 
- { (q - A) (2* - 2q + 2A + f ) Cj (S 2 r x r y ) } f q \ 

+ {{q-\){q-\-l)/(S 2 rl)}f^_ 2 . 



where r x , r y and oj are constant. 

Due to the spherical symmetry, the energy level with the given total 
angular momentum £ is (2£ + l)-degeneracy (normal degeneracy). Further- 
more, since A < q < £, there are (^+1) sets of solutions with the parity (— 1) 
and £ sets of solutions with the parity (— l) e+1 . This conclusion coincides 
with that by Wigner (see P. 218 in [17]). When t = we have the constant 
solution with zero energy and even parity. When I = 1, we have one set of 
solutions ^e m i with the even parity and two sets of solutions ^e m o with the 
odd parity: 



It is similar to obtain the solutions with the higher orbital angular momen- 
tum £. The partners of the solutions with the smaller eigenvalues of L3 can 
be calculated from them by the lowering operator L_. 

4. CONCLUSION 

In summary, we have reduced the three-body Schrodinger equation for any 
given total orbital angular momentum and parity to a system (12) of the 
coupled partial differential equations with respect only to three variables, 
describing the internal degrees of freedom in a three-body problem. This 
equation system is a good starting point for solving the quantum three-body 
problems. As an example, we obtain the exact solution to the Schrodinger 
equation for the rigid body. 
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(22) 



E$ =h 2 /(2Mr 2 y sm 2 uj) , 




) 



(2/1 +m), 



(23) 
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